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Nuclear collisions recreate conditions in the universe microseconds after the 
Big Bang. Only a very small fraction of the emitted fragments are light nu- 
clei, but these states are of fundamental interest. We report the observation 
of antihypertritons - composed of an antiproton, antineutron, and antilambda 
hyperon - produced by colliding gold nuclei at high energy. Our analysis yields 
70 ± 17 antihypertritons (|H) and 157 ± 30 hypertritons (aH). The measured 
yields of (|H) and ^He (^He) are similar, suggesting an equilibrium in co- 
ordinate and momentum space populations of up, down, and strange quarks 
and antiquarks, unlike the pattern observed at lower collision energies. The 
production and properties of antinuclei, and nuclei containing strange quarks, 
have implications spanning nuclear/particle physics, astrophysics, and cosmol- 
ogy. 

Nuclei are abundant in the universe, but antinuclei that are heavier than the antiproton have 
been observed only as products of interactions at particle accelerators (7, 2). Collisions of 
heavy nuclei at the Relativistic Heavy-Ion Collider (RHIC) at Brookhaven National Laboratory 
(BNL) briefly produce hot and dense matter that has been interpreted as a quark gluon plasma 
(QGP) (3, 4) with an energy density similar to that of the universe a few microseconds after 
the Big Bang. This plasma contains roughly equal numbers of quarks and antiquarks. As a 
result of the high energy density of the QGP phase, many strange-antistrange {ss) quark pairs 
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are liberated from the quantum vacuum. The plasma cools and transitions into a hadron gas, 
producing nucleons, hyperons, mesons, and their antiparticles. 

Nucleons (protons and neutrons) contain only up and down valence quarks, while hyperons 
(A, S, S, f^) contain at least one strange quark in its 3 -quark valence set. A hypernucleus is a 
nucleus that contains at least one hyperon in addition to nucleons. All hyperons are unstable, 
even when bound in nuclei. The lightest bound hypernucleus is the hypertriton (aH), which 
consists of a A hyperon, a proton, and a neutron. The first observation of any hypernucleus 
was made in 1952 using a nuclear emulsion cosmic ray detector (5). Here, we present the 
observation of an antimatter hypernucleus. 

Production of antinuclei: Models of heavy-ion collisions have had good success in ex- 
plaining the production of nuclei by assuming that a statistical coalescence mechanism is in 
effect during the late stage of the collision evolution (4, 6). Antinuclei can be produced through 
the same coalescence mechanism, and are predicted to be present in cosmic rays. An observed 
high yield could be interpreted as an indirect signature of new physics, such as Dark Mat- 
ter (7, 8). Heavy-ion collisions at RHIC provide an opportunity for the discovery and study of 
many antinuclei and antihypernuclei. 

The ability to produce antihypernuclei allows the study of all populated regions in the 3- 
dimensional chart of the nuclides. The conventional 2-dimensional chart of the nuclides orga- 
nizes nuclear isotopes in the (N, Z) plane, where is the number of neutrons and the Z is 
the number of protons in the nucleus. This chart can be extended to the negative sector in the 
(A^, Z) plane by including antimatter nuclei. Hypernuclei bring a third dimension into play, 
based on the strangeness quantum number of the nucleus. The present study probes the terri- 
tory of antinuclei with non-zero strangeness (Fig. 1), where proposed ideas (9-12) related to 
the structure of nuclear matter can be explored. 

Hypernuclei-Formation and observation: The hyperon-nucleon (YN) interaction, respon- 
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Figure 1 : A chart of the nuchdes showing the extension into the strangeness sector. Normal 
nuclei lie in the (A^, Z) plane. Antinuclei lie in the negative sector of this plane. Normal 
hypernuclei lie in the positive (A^, Z) quadrant above the plane. The antihypertriton |H reported 
here extends this chart into the strangeness octant below the antimatter region in the (A^, Z) 
plane. 

sible in part for the binding of hypernuclei, is of fundamental interest in nuclear physics and 
nuclear astrophysics. For example, the YN interaction plays an important role in attempts to 
understand the structure of neutron stars. Depending on the strength of the YN interaction, the 
collapsed stellar core could be composed of hyperons, of strange quark matter, or of a kaon 
condensate {13). While the hyperons or strange particles inside a dense neutron star would not 
decay because of local energy constraints, free hypernuclei decay into ordinary nuclei with typ- 
ical lifetimes of a few hundred picoseconds, which is still thirteen orders of magnitude longer 
than the lifetimes of the shortest-lived particles. The lifetime of a hypernucleus depends on the 
strength of the YN interaction (74, 75). Therefore, a precise determination of the lifetime of 
hypernuclei provides direct information on the YN interaction strength (75, 16). 

The experiment was carried out by the STAR collaboration (77) at the RHIC facility. The 
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main detector of the STAR experiment is a gas-filled cylindrical Time Projection Chamber 
(TPC), with an inner radius of 50 cm, an outer radius of 200 cm, and a length of 420 cm along 
the beam line (18). The TPC is a device for imaging, in three dimensions, the ionization left 
along the path of charged particles. It resolves over 50 million pixels within its active volume. 
The present analysis is based on interactions produced by colliding two Au beams at an energy 
of 200 GeV per nucleon-nucleon collision in the center-of-mass system. Approximately 89 
million collision events were collected using a trigger designed to accept, as far as possible, all 
impact parameters (minimum-bias event), and an additional 22 million events were collected 
using a trigger that preferentially selects near-zero impact parameter (or "head-on") collisions. 
The accepted collisions are required to occur within 30 cm of the center of the TPC along 
the beam line. Charged particle tracks traversing the TPC are reconstructed in an acceptance 
that is uniform in azimuthal angle. The precise coverage in terms of polar angle is somewhat 
complicated (18), but roughly speaking, charged tracks emerging at angles with respect to the 
beam axis in the range of 45° < 6 < 135° are reconstructed. 

Fig. 2 depicts a typical Au+Au collision reconstructed in the STAR TPC. The tracks are 
curved by a uniform magnetic field of 0.5 Tesla parallel to the beam line. The event of interest 
here includes a |H candidate created at the primary collision vertex near the center of the TPC. 
The |H travels a few centimeters before it decays. One of the possible decay channels is 
^He+7r+, which occurs with a branching ratio of 25% assuming that this branching fraction is 
the same as that for (75). The two daughter particles then traverse the TPC along with the 
hundreds of other charged particles produced in the primary Au+Au collision. The trajectories 
of the daughter particles are reconstructed from the ionization trails they leave in the TPC gas 
volume (shown in Fig. 2 as thick red and blue lines for ^He and 7r+, respectively). The energy 
loss by these particles to ionization in the TPC, {dE/dx), depends on the particle velocity 
and charge. Particle identification is achieved by correlating the {dE/dx) values for charged 
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Figure 2: A typical event in the STAR detector that includes the production and decay of a 
|H candidate. In (A), the beam axis is normal to the page, and in (B), the beam axis is hori- 
zontal. The dashed black line is the trajectory of the |H candidate, which cannot be directly 
measured. The heavy red and blue lines are the trajectories of the ^He and 7r+ decay daughters, 
respectively, which are directly measured. 

particles in the TPC with their measured magnetic rigidity, which is proportional to the inverse 
of the curvature of the trajectory in the magnetic field. With both daughter candidates directly 
identified, one can trace back along the two helical trajectories to the secondary decay point, 
and thereby reconstruct the location of the decay vertex as well as the parent momentum vector. 

Particle identification: Fig. 3 presents results from the antihypertriton analysis outlined 
above, along with results from applying the same analysis to measure the normal matter hy- 
pertritons in the same dataset — only the sign of the curvature of the decay products is re- 
versed. Fig. 3C shows (dE/dx) for negative tracks as a function of the magnetic rigidity; 
the different bands result from the different particle species. The measured (dE/dx) of the 
particles is compared to the expected value from the Bichsel function (79), which is an ex- 
tension of the usual Bethe Bloch formulas for energy loss. A new variable, z, is defined as 
z = ln{{dE/dx) / (dE/dx) ^), where (dE/dx) b is the expected value of (dE/dx) for the given 
particle species and momentum. The measured z(^He) distributions for ^He and ^He tracks 
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(Fig. 3D), includes 5810 ^He and 2168 ^He candidates with |z(^He)| < 0.2, and represents the 
largest sample of ^He antinuclei that has been collected to date. The first few ^He candidates 
were observed at the Serpukhov accelerator laboratory (20), followed by confirmation from the 
European Organization for Nuclear Research (CERN) (27). In 2001, a relatively large ^He sam- 
ple was reported by the STAR collaboration (22). The ^He and ^He samples in the present 
analysis are so cleanly identified that misidentification from other weak decays is negligible. 
However, due to the (dE/dx) overlap between and ^He at low momenta, it is only possible 
to identify the ^He nuclei at relatively high momenta (i.e. above ^2 GeV/c). The daughter 
pions from |H decays usually have momenta ^ 0.3 GeV/c, and can be cleanly identified (23). 

Topological reconstruction: A set of topological cuts is invoked in order to identify and 
reconstruct the secondary decay vertex positions with a high signal-to-background ratio. These 
cuts involve the distance at the decay vertex between the tracks for the ^He and 7r+ (< 1 cm), the 
distance of closest approach (DCA) between the |H candidate and the event primary vertex (<1 
cm), the decay length of the |H candidate vertex from the event primary vertex (>2.4 cm), and 
the DCA between the tt track and the event primary vertex (>0.8 cm). The cuts are optimized 
based on full detector response simulations (24). Several different cut criteria are also applied 
to cross-check the results and to estimate the systematic errors. The signal is always present, 
and the difference in the total yields using different cuts are found to be less than 15%. The 
total systematic error in the present analysis is 15%. 

The parent candidate invariant mass is calculated based on the momenta of the daughter 
candidates at the decay vertex. The results are shown as the open circles in Fig. 3A for the 
hypertriton: -^^ He + tt", and in Fig. 3B for the antihypertriton: |H ^He+7r+. There 
remains an appreciable combinatorial background in this analysis, which must be described 
and subtracted. A track rotation method is used to reproduce this background. This approach 
involves the azimuthal rotation of the daughter ^He (^He) track candidates by 180 degrees with 
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Figure 3: (A, B) show the invariant mass distribution of the daughter ^He + tt. The open circles 
represent the signal candidate distributions, while the solid black lines are background distri- 
butions. The blue dashed lines are signal (Gaussian) plus background (double exponential) 
combined fit (see the text for details). A (B) shows the (|H) candidate distributions. (C) 
shows (dE/dx) versus rigidity (momentum/ 1 nuclear charge units |) for negative tracks. Also 
plotted are the expected values for ^He and tt tracks. (D) and (C) demonstrate that the ^He and 
^He tracks (|z(^He)| < 0.2) are identified essentially without background. 

respect to the event primary vertex. In this way, the event is not changed statistically, but all of 
the secondary decay topologies are destroyed because one of the daughter tracks is rotated away. 
This provides an accurate description of the combinatorial background. The resulting rotated 
invariant mass distribution is consistent with the background distribution, as shown by the solid 
histograms (Fig. 3A,B). The rotated background distribution is fit with a double exponential 
function: f{x) oc exp(— ^) — exp(— ^), where x = m — m(^He) — m(7r), andpi, p2 are fit 
parameters. Finally, the counts in the signal are calculated after subtraction of this fit function 
derived from the rotated background. In total, 157 ± 30 and 70 ± 17 |H candidates are thus 
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observed. The quoted errors are statistical. 

Production and properties: We can use the measured yield to estimate the expected 
yield of the |H, assuming symmetry between matter and antimatter, in the following manner: 
|H = ^Hx^He / ^He= 59 ± 11. This indicates a 5.2a projection of the number of |H that is 
expected in the same data set where ^H, ^He and ^He are detected. An additional check involves 
fitting the ^He +7r invariant mass distribution with the combination of a Gaussian "signal" term 
plus the double-exponential background function (the blue-dashed lines in Fig. 3A,B). The 
resulting mean values and widths of the invariant mass distributions are consistent with the 
results from the full detector response simulations. Our best fit values (from ^ minimization) 
arem(iH)= 2.989±0.001±0.002GeV/c2 and m(|H) = 2.991 ±0.001 ±0.002 GeV/c^. These 
values are consistent with each other within the current statistical and systematic errors, and are 
consistent with the best value from the literature, i.e., m(^H) = 2.99131 ±0.00005 GeV/c^ (76). 
Our systematic error of 2 MeV/c^ arises from well-understood instrumental effects that cause 
small deviations from ideal helical ionization tracks in the TPC. 

Lifetimes: The direct reconstruction of the secondary decay vertex in this data allows mea- 
surement of the lifetime, r, via the equation Nit) = A^(0)e"^/^, t = l/{(5^c), f3^c = p/m, 
where / is the measured decay distance, p is the particle momentum, m is the particle mass, and 
c is the speed of light. For better statistics in our fit, the and |H samples are combined, as 
the matter-antimatter symmetry requires their lifetimes to be equal. Separate measurements of 
the lifetimes for the two samples show no difference within errors. The signal is then plotted in 
three bins in 1/I3j. The yield in each bin is corrected for the experimental tracking efficiency 
and acceptance. The total reconstruction efficiency of the |H and is on the order of 10%, 
considering all sources of loss and the analysis cuts. The three points are then fit with the expo- 
nential function to extract the parameter cr, and the best-fit result is displayed as the solid line 
in Fig. 4A. To arrive at the optimum fit, a analysis was performed (see the inset to Fig. 4A). 
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Figure 4: (A) The (solid squares) and A (open circles) yield distributions versus cr. The 
solid lines represent the cr fits. The inset depicts the distribution of the best cr fit. (B) 
World data for lifetime measurements. The data points are from Refs. (25-30). The theo- 
retical calculations are from Refs. (14, 15). The error bars represent the statistical uncertainties 
only. 

The CT parameter that is observed in this analysis is cr = 5.5 ±^ 4 ±0.8 cm, which corresponds 
to a lifetime r of 182 ±45 ±27 ps. As an additional cross-check, the A hyperon lifetime is 
extracted from the same data set using the same approach, for the A ^ p ± tt" decay channel. 
The result obtained is r = 267 ± 5{stat) ps, which is consistent with r = 263 ± 2 ps compiled 
by the Particle Data Group (19). 

The lifetime measurements to date (25-31) are not sufficiently accurate to distinguish 
between models, as depicted by Fig. 4B. The present measurement is consistent with a calcula- 
tion using a phenomenological wave function (74), and is also consistent with a more recent 
three-body calculation (75) using a more modern description of the baryon-baryon force. The 
present result is also comparable to the lifetime of free A particles within the uncertainties, and 
is statistically competitive with the earlier experimental measurements. 

Coalescence calculations: The coalescence model makes specific predictions about the ra- 
tios of particle yields. These predictions can be checked for a variety of particle species. To 
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Table 1 : Particle ratios from Au+Au collisions at 200 GeV. 



Particle type 



Ratio 




0.49 ±0.18 ±0.07 



^He/^He 



0.45 ± 0.02 ± 0.04 



iH/^He 



0.89 ±0.28 ±0.13 



jH/^He 



0.82 ±0.16 ±0.12 



determine the invariant particle yields of |H and ^H, corrections for detector acceptance and 
inefficiency are applied. The |H and yields are measured in three different transverse mo- 
mentum (pt) bins within the analyzed transverse momentum region of 2 < < 6 GeV/c and 
then extrapolated to the unmeasured regions (pt < 2 GeV/c and Pt > Q GeV/c). This extrapo- 
lation assumes that both |H and have the same spectral shape as the high-statistics ^He and 
^He samples from the same data set (see Table 1). If the |H and are formed by coalescence 
of (A ± p + n) and (A + p + n), then the production ratio of |H to should be proportional to 

X ^ X f). The latter value can be extracted from spectra already measured by STAR, and the 
value obtained is 0.45 ± 0.08 ± 0.10 (23, 24), The measured |H / and ^He / ^He ratios are 
consistent with the interpretation that the |H and are formed by coalescence of (A + p + n) 
and (A + p ± n), respectively. 

Discussion: As the coalescence process for the formation of (anti)hypernuclei requires that 
(anti)nucleons and (anti)hyperons be in proximity in phase space (i.e., in coordinate and mo- 
mentum space), (anti)hypernucleus production is sensitive to the correlations in phase-space 
distributions of nucleons and hyperons (6). An earlier two-particle correlation measurement 
published by STAR implies a strong phase-space correlation between protons and A hyper- 
ons {32). Equilibration among the strange quark flavors and light quark flavors is one of the 
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proposed signatures of QGP formation (33), which would result in high (anti)hypernucleus 
yields. In addition, recent theoretical studies motivate a search for the onset of QGP by studying 
the evolution of the baryon - strangeness correlation as a function of collision energy (34-36). 
The yields provide a natural and sensitive tool to extract this correlation (6, 37), as they can 
be compared to the yields of ^He and ^H, which have the same atomic mass number. Besides 
4m + 4d valence quarks, the valence quark content of these species includes one additional u, d 
and s quark for ^He, and ^H, respectively. Recent nuclear transport model calculations (37) 
support the expectation that the strangeness population factor,5'3 =|H/(^He xA/p), can be used 
as a tool to distinguish the QGP from a purely hadronic phase. 
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Figure 5: Particle ratios as a function of center-of-mass energy per nucleon-nucleon collision. 
The data at lower energies are from Refs. (38^0). The ^H/(^He x A/p) ratio is corrected for 
the spin degeneracy factor (38). The error bars represent statistical uncertainties only. 



Fig. 5 depicts various particle ratios as a function of the collision energy. The ^He/^H ratio 
at a center-of-mass energy of 5 GeV obtained at the Alternating Gradient Synchrotron (AGS) 
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at BNL is much closer to unity than the ratio ^H/^He at the same energy. The values of 6^3 are 
about 1/3 at AGS energies, and near unity at RHIC energies, although with large uncertainties. 
The AGS value is further constrained to be relatively low by the measured upper limit on the 
^H/^He ratio (38), indicating that the phase space population for strangeness is very similar to 
that for the light quarks in high-energy heavy-ion collisions at RHIC, in contrast to the situation 
at AGS. 

Individual relativistic heavy-ion collisions produce abundant hyperons containing one (A, S), 
two (S) or three (Q) strange (anti)quarks. The coalescence mechanism for hypernucleus pro- 
duction in these collisions thus provides a source for other exotic hypernucleus searches. This 
should allow an extension of the 3-D chart of the nuclides (Fig. 1) further into the antimat- 
ter sectors. Future RHIC running will provide increased statistics, allowing detailed studies 
of masses and lifetimes, as well as stringent tests of production rates compared to predictions 
based on coalescence models. 

Concluding remark: Evidence for the observation of an antihypernucleus, the |H, with 
a statistical significance of 4. la has been presented; consistency checks and constraints from 
a analysis in the same event sample, with 5.2a significance, support this conclusion. The 
lifetime is observed to be r = 182 ±45 ±27 ps, which is comparable to that of the free A hyperon 
within current uncertainties. The |H (^H) to ^He (^He) ratio is close to unity and is significantly 
larger than that measured at lower beam energies, indicating that the strangeness phase space 
population is similar to that of light quarks. An order-of-magnitude larger sample of similar 
collisions is scheduled to be recorded in the near future. The antihypernucleus observation 
demonstrates that RHIC is an ideal facility for producing exotic hypernuclei and antinuclei. 
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